Mechanism for flux guidance by micrometric antidot arrays in superconducting films 



J. I. Vestgarden, 1 V. V. Yurchenko, 1 R. Wordenweber, 2 and T. H. Johansen 1 ' 3,4 

1 Department of Physics, University of Oslo, P. 0. Box 1048 Blindern, 0316 Oslo, Norway 
2 Peter Grilnberg Institute and JAR A- Fundamentals of Future Information Technology, 
Forschungszentrum Jillich, 52425 Jiilich, Germany 
3 Institute for Superconducting and Electronic Materials, 
University of Wollongong, Northfields Avenue, Wollongong, NSW 2522, Australia 
4 Centre for Advanced Study at The Norwegian Academy of Science and Letters, Drammensveien 78, 0271 Oslo, Norway 

A study of magnetic flux penetration in a superconducting film patterned with arrays of micron 
sized antidots (microholes) is reported. Magneto-optical imaging (MOI) of a YBa2Cu3 x film 
shaped as a long strip with perpendicular antidot arrays revealed both strong guidance of flux, and 
at the same time large perturbations of the overall flux penetration and flow of current. These 
results are compared with a numerical flux creep simulation of a thin superconductor with the 
same antidot pattern. To perform calculations on such a complex geometry, an efficient numerical 
scheme for handling the boundary conditions of the antidots and the nonlocal electrodynamics was 
developed. The simulations reproduce essentially all features of the MOI results. In addition, the 
numerical results give insight into all other key quantities, e.g., the electrical field, which becomes 
extremely large in the narrow channels connecting the antidots. 

PACS numbers: 74.25.Ha, 74.25.Op 



I. INTRODUCTION 

The motion of magnetic flux in type-II supercon- 
ducting films can to a large extent be controlled by 
introduction of artificial micro- and nano- structures, 
such as antidots (holes) POD magnetic dotsf^ thick- 
ness modulations grain boundaries slits or 
magnet ic do main walls in superconductor/ferromagnet 
hybrids (HEU Such structures are key building blocks 
for successful realization of fluxonics devices like vor- 
tex ratchets, pumps and lenses etcP^^It is known that 
when antidots are sufficiently small they can become pin- 
ning sites for the vortices This type of pinning is usu- 
ally noticable only close to the critical temperature T c , 
and can be observed, e.g., as pronounced matching be- 
tween t he vo rtex density and the underlying lattice of 
antidots! 17 ^ 18 " It was also demonstrated that certain pat- 
terns can be used to reduce noise due to vortex motion 
in SQUIDsp^l and quite recently in superconducting mi- 
crowave resonators P3 

Of technological as well as fundamental interest is also 
realizations of flux guidance, i.e., how to achieve directed 
and controlled motion of the magnetic flux. When the 
pinning by the antidots is strong compared to the intrin- 
sic pinning of the material, guidance is effectuated by the 
more mobile interstitial vorticesP^ Conversely, when the 
intrinsic pinning is strong, flux moves most easily inside 
the holes, and arrays of antidots should enhance the flux 
penetration.^ It has bee n sho wn that by introducing ar- 
rays of dotP^ or antidot $£&M VO rtices can be guided 
away from the direction given by the Lorentz force of an 
applied current. Similarly, a periodic arrangement of an- 
tidots can cause effective flux drainage of a sample in the 
descending branch of a magnetic field ramp It was also 
found that when the antidot lattice breaks the symmetry 
of the overall sample shape, or when the antidots have 



nontrivial shap es, the critical current density can become 
anisotropic 3 -^ and the flux motion can be enhanced in 
unexpected directions.^ 

Whereas superconducting films containing complex ar- 
rangements of antidots can today be readily produced 
using, e.g., optical lithography, the theoretical modelling 
of the local and global features of the flux dynamics 
in such systems is challenging. In essence, this is due 
to the nonlocal nature of the electrodynamics in two- 
dimensional samples subjected to a perpendicular mag- 
netic field. In this case the critical-state is strongly mod- 
ified since at any applied field there will be currents flow- 
ing over the whole sample area, including regions in the 
flux-free Meissner-state! 28 " 30 * Also numerical simulations 
of the flux dynamics are a lot more computationally de- 
manding in films compared to bulk.^^ Additional com- 
plications appear in a film with antidots, as the shield- 
ing currents then meet constrictions, and the flow must 
adapt to the available and often narrow bridges of super- 
conducting material between the antidots. Hence, the 
current density quickly rises to the critical value, causing 
major rearrangements of both the current flow and the 
flux distribution.^ 

Magneto-optical imaging (MOI) is a unique experi- 
mental tool for observing the nontrivial redistributions 
created by patterning of superconducting films. The 
technique allows direct observation of the flux density 
over length scales ranging from the entire sample size 
and down to the size of individual antidots. In spite 
of large previous efforts, surprisingly few investigations 
were carried out on films with simple antidot patterns. 
In this work we present results from MOI experiments to- 
gether with theoretical modelling of the flux and current 
behavior in a thin superconducting strip containing only 
a few linear arrays of antidots. The experiments were 
performed using a film of YE^CusO^ (YBCO) cooled 



far below T c to be in the regime where flux pinning by 
the antidots is negligible. 

The paper is organized as follows. In Sec. II we outline 
the current flow modification caused by an antidot array 
as expected from the critical state model. Section III 
presents results from our MOI experiments. Section IV 
describes our method for numerical simulations of the 
electrodynamics, and the results are presented and dis- 
cussed in Sec. V. Section VI gives a summary. 



II. CRITICAL CURRENT FLOW LINES 

Before presenting the experimental results we con- 
sider the predictions of the Bean critical-state model for 
a type-II superconducting strip with a linear array of 
equally spaced antidots crossing from one side to the 
other, see Fig. [I] The antidots have radius r and center- 
to-center distance a. An external magnetic field, H a , is 
applied perpendicular to the strip plane, and we assume 
the superconductor initially contains no flux and current. 

As the field is gradually increased magnetic flux will 
penetrate from the edge and shielding currents will be- 
gin to flow. In areas containing flux, the critical state is 
formed and the current density has the critical magni- 
tude j c , while in the unpenetrated Meissner-state region 
the currents have a subcritical density. To draw the cur- 
rent stream lines in the critical-state region we use the 
simple rule following from the Bean critical-state model, 
namely that the magnitude of the current density is con- 
stant, and thus represented by a set of equidistant stream 
lines with spacing inversely proportional to j c . The con- 
struction starts by drawing continuous lines that follow 
the external perimeter of the sample, and when reach- 
ing a hole the current must flow around it. This scheme 
results in the flow pattern shown in the Fig. [T] corre- 
sponding to a partly penetrated state. The closure of 
each streamline takes place in parts of the sample not 
included in the figure. 

The figure shows that the presence of the antidot array 
deforms the critical current in a large region with a near 
rhombic shape. A line denoted d-linei marks where the 
direction of the current flow suddenly begins to deviate 
from being parallel to the edge. Using the principle of 
current flow continuity one finds that the line makes an 
angle a with the edge, given by 



cos a 



(i) 



Evidently, such a line exists on both sides of the array, 
and the pair defines one half of the rhombic area. 

In addition to the large perturbation of the current 
flow, the construction also results in a fine structure in 
the flow pattern. Due to the shape of the antidots the 
stream lines inside the rhombic area consist of a sequence 
of circular arcs. As seen from the figure, each antidot 
creates its own pair of lines along a direction defined by 
the cusped joints of two arcs. These lines, denoted d- 
line2, become straight a short distance from the array, 




d-line. 



FIG. 1. (Color online) Construction of the critical current 
streamlines near a linear array of equidistant circular antidots. 
The pattern follows from the Bean critical-state model for a 
sample with a straight edge. Note that streamlines of the 
sub-critical currents flowing in the Meissner-state part of the 
strip are not included. 



and one finds that the angle, (3, they make with the strip 
edge is the same for all of them, and given by 



cot/3 



, yV(a-r) 
a-2r 



(2) 



Related to this construction is a previous analysis of 
the current flow in the case of having a weak link across 
the superconducting strip P It was assumed that inside 
the weak link the critical current density is uniformly 
reduced by a factor jb/j c - Also that construction re- 
sulted in a rhombic area where the current has changed 
direction, and the angle a was found to be given by 
cos 2a = —(jb/jc)- Indeed, this is equivalent to Eq. (fi]) 
with the transparency of the antidot array, (a — 2r)/a, 
replaced by the current ratio in the weak link case. 



III. EXPERIMENT 

A 150 nm thick YBCO film was produced by mag- 
netron sputtering on a r-cut sapphire substrate. The 
sample was shaped using optical lithography and ion 
beam etching into a long strip 0.5 mm in width, and with 
an arrangement of antidots as shown in Fig. [2^i. The an- 
tidots having radius r = 2 /im form a linear array with 
period a = 10 /im. Along the center line there are two 
larger antidots of radius 20 /im, connected with an array 
of 6 antidots, also of radius r = 2 /im. The arrangement 
was motivated by a possibility to study flux guidance 
along linear antidot arrays both perpendicular and par- 
allel to the strip edges. The large holes are included to 
serve as reservoirs for the incoming flux. 
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Magneto-optical imaging of the sample was performed 
using a bismuth substituted ferrite garnet film with in- 
plane magnetization as Faraday rotating sensor Primages 
of the flux distribution were recorded through a polarized 
light microscope using crossed polarizers. In this way the 
image brightness represents the magnitude of the flux 
density. For details of the setup, see Ref. |35] 



x 




FIG. 2. (Color online) Sketch of the antidot arrangement 
in the superconducting sample, which has an overall shape as 
a long thin strip in the xy-plane (a). The flux distribution, 
B z , in the YBCO film seen by MOI and recorded at applied 
fields of 2.4 mT, 5.7 mT and 10 mT, respectively (b-d). 




FIG. 3. (Color online) Details of Fig. [2] showing the d- 
line pattern from a row of antidots. The straight d-linei and 
d-line2 are drawn as a guide to the eye with angles corre- 
sponding to Eqs. |l]) and Q calculated for the geometrical 
characteristics of the antidot array in the YBCO film. 



Shown in Fig. [2}>d is the flux distribution, B Z1 in the 
patterned part of the sample as the applied field is slowly 
ramped up to iioH a = 10 mT after an initial zero-field- 
cooling to 50 K. As typical for thin films, the rim of the 
strip appears bright, thus showing piling up of the field 
that is expelled by the superconductor. Already at a field 
of 2.4 mT (panel b) the linear arrays perpendicular to 
the edges are clearly visible, giving direct evidence that 
under these conditions the antidots are guiding the flux 
rather than pinning the vortices. When the field becomes 
5.7 mT (panel c) guided flux reaches the large hole in the 
center of the strip, and interestingly, the flux continues 
its motion out of the hole being guided by the antidot 
array oriented parallel to the edges. In panel d, one sees 
that at 10 mT the parallel guidance has filled also the 
second large hole with flux well before the strip is fully 
penetrated. 

Comparing the MOI results with the current flow pat- 
tern outlined in Fig. [I] note in the image of Fig. [2]i that 
four straight dark lines form a diamond shape where the 
two transverse antidot arrays make up the long diago- 
nal. These dark lines correspond to discontinuity lines of 
the type marked as d-linei in the critical-state construc- 
tion. From Eq. and the dimensional characteristics of 
the present antidot array the angle a should be given by 
cos a = 1/a/5, or a = 63°. Shown in Fig. [3] is a close-up 
view of B z near the antidot array superimposed with a 
pair of dotted lines tilted by a = 63°, demonstrating an 
excellent quantitative agreement. Also the fact that the 
diamond- feature shows up as dark, i.e., has a very low 
flux density, follows naturally from the construction in 
Fig. [l] since the sharp clock- wise turning of the current 
near the d-linei provides an additional local shielding. 

The fine structure of B z in the region around the an- 
tidot arrays is only faintly visible. This is not surprising 
considering the d-line2 in Fig. [I] where the turning of the 
currents at the cusps is gradually reduced away from the 
antidots. Nevertheless, one can clearly see in Fig.[3]traces 
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of a fishbone-pattern. The angle (3 should according to 
Eq. (|| be /? = 37°, and Fig. [3] includes as a guide to the 
eye a pair of dotted lines having this angle. One finds a 
quite nice agreement with the streaky features visible in 
the magneto-optical image on both sides of the antidot 
array. 

Note that also the observed high flux density along 
the array of antidots can be readily understood from the 
current flow construction in Fig. [T] As the current flows 
past the constricted region the stream lines make a very 
sharp turn in the counter-clockwise direction. This cur- 
vature enhances flux density of the same polarity as the 
one entering from the edge, thus creating an effective flux 
guidance well beyond the overall flux penetration front 
in the strip. 

Although many features of the MOI results can be un- 
derstood from the critical-state considerations above, this 
picture is far from complete. In particular, the entire dis- 
tribution of currents flowing in the Meissner-state part of 
the film was neglected. Moreover, the electric field is not 
considered in the simple analysis, and is also not available 
from the experiment. 

To complete the analysis of the flux guidance, we have 
developed an efficient numerical scheme allowing to carry 
out simulations of the electrodynamic response of a thin 
superconductor with antidots. Below is a description of 
the scheme followed by a report of the numerical results 
for a film patterned just like the present YBCO sample. 



IV. SIMULATION SCHEME 

The numerical scheme assumes that the superconduc- 
tor is thin, i.e., it has a thickness d much less than any 
lateral dimension of the sample. The external field, H ai 
is applied in the perpendicular z-direction, and the in- 
duced currents will be quantified by the sheet current, J. 
When J approaches the critical magnitude J c = dj c , the 
depinning transition is sharp, and gives rise to a highly 
nonlinear material characteristics conventionally approx- 
imated by a power law, 



E = pj/d, p = p (J/JcY 



(3) 



where E is the electric field and p is the resistivity with 
po being a characteristic value. The creep exponent n 
is usually large, typically in the range 20 < n < 60 for 
YBCOP^ The Bean model corresponds to the limit 
n — » oo. 

Rather than working directly with the sheet current, 
a more convenient quantity is the local magnetization 
g = g(r,i), defined by 31 



dg_ 
dy 



— Jxi 



(4) 



where r = (x,y). The g incorporates current conserva- 
tion since it by definition gives V • J = 0. The function 



g is extended to the whole space by setting g = out- 
side the sample. Inside the antidots g is uniform, but not 
necessarily zero. 

Another basic equation is the Biot-Savart law, which 
can be expressed as 



B z /po = H a + Qg, 
with the operator Q given b}E^ 



Q9(r) = T~ x 



(5) 



(6) 



where T is the 2D spatial Fourier transform, and k = |k|. 
The inverse relation is 



"VW=J- 1 



(7) 



where cp is an auxiliary function. 

By taking the time derivative of Eq. (J5|, we get 
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B z /po — H a 



(8) 



which is solved by discrete integration forward in time. 
This is possible since B z (r,t) can be calculated from 
#(r,t), as described below. 

To solve Eq. ([8| the space is discretized in such way 
that Q and its inverse can be implemented using Fast 
Fourier Transforms. We let the superconductor occupy 
the space \x\ < w and y < \w\, and by using periodic 
boundary conditions along the sample has the shape 
of a long strip in the ^/-direction. In the x-direction we let 
the strip be surrounded by empty space so that the total 
area included in the calculations is \x\ < 1.5k;, y < \w\. 

Hence, the x?/-plane consists of three different parts: 
the superconductor, the area outside the strip, and the 
area inside the antidots. In order to solve Eq. ([8|, we 
must find B z in all three regions, which requires different 
algorithms. 

Starting with the superconductor itself, it obeys the 
material law, Eq. ([3|, which when combined with Fara- 
day's law, B z = —(V x E) 2 , gives 



B Z =V -(pVg)/d. 



(9) 



From g(r,t) the gradient is readily calculated, and since 
the result allows finding J(r, t) from Eq. Q, also p(r,t) 
is determined from Eq. The task then is to find B z 
in the non-superconducting parts, so that g = outside 
the strip and g is uniform in the antidots. This cannot 
be calculated efficiently using direct methods due to the 
nonlocal relation between B z and g. Instead we use an 
iterative procedure. 

For all iteration steps, i = l...s, B z ^ is fixed inside the 
superconductor by Eq. At i = 1, an initial guess is 

made for B^ outside the sample and inside the antidots, 
and (jM is calculated from Eq. In general, this g^ 
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FIG. 4. (Color online) Numerical calculation of the distribution of B z in a superconducting thin strip after an applied magnetic 
field was increased from zero to H a = 1. The strip contains a pattern of antidots essentially identical to the sample studied 
experimentally. As in MOI the image brightness represents the magnitude of B z as indicated in the color bar insert. 



does not vanish outside the strip and an improvement is 
obtained by 



(10) 



The projection operator O is unity outside the strip and 
zero everywhere else. Also, the output of the operation 
should be shifted to satisfy j cPrOg = 0. The constant 
is determined by requiring flux conservation, 



d 2 r [B$ 



Correspondingly, B z 1 ^ 1 ^ 



(11) 



in each of the antidots is also 
found by Eq. fllOl, but where the projection operator O 



now is unity in the antidot and zero everywhere else, and 
with calculated using Faraday's law, 



i as. 



dl-E 



antidot boundary 



d 2 rB z . (12) 



antidot area 



Thus, at each iteration (i), we run through all antidots 

and the outside area, and calculate B z l+1 \ The proce- 
dure is repeated until after i = s iterations becomes 
sufficiently uniform both outside the strip and within the 
antidots. Then, g^ is inserted in Eq. ([8]), which brings 



us to the next time step, where the whole iterative pro- 
cedure starts anew. 

Several comments can be made regarding the imple- 
mentation of the simulation scheme. First, the algorithm 
for finding g scales as 0(N\ogN) with the total number 
of discrete grid points N. This enables simulations of 
large grids, say N > 10 6 . Second, the quality and per- 
formance is improved by a good initial value, such as 



B^\t) = B z S) (t - At). Third, for small antidots there 
is a large perfomance gain by replacing the full operator 
QO in Eq. (10) by a local operator that only runs over 



- rW 



the antidot. Fourth, the results can be made more robust 
by enforcing J = exactly by hand after each iteration 
step to prevent accumulations of small unphysical cur- 
rents inside the antidots with time. Finally, a robust way 
to satisfy Eq. ( 12 ) is to first calculate B z from Eq. ([9| in 
the antidot area using the material law as if the antidot 
was not there, and then choose so that the total flux 
in the antidot is the same. 

At all times, the simulation scheme provides direct ac- 
cess to the distributions of g, B Z1 J, everywhere in the 
plane z = 0, as well as access to E inside the supercon- 
ductor, through Eq. (J3|. Note that E is not calculated 
inside the antidots or outside the sample. 



V. SIMULATION RESULTS 

Numerical simulations were carried out for a supercon- 
ducting strip of width 2w where between the two edges 
there is a linear array of 40 antidots of radius r — 0.01k; 
and center-to-center separation a = 4.4r. Two larger 
holes of radius 0.08iu are located on the strip center 
line, and with another array of 6 antidots in between. 
The area of calculation, 3w x 2w, was discretized on a 
1536 x 1024 equidistant grid. The creep exponent was 
n = 29. The simulation was carried out in dimension- 
less units based on J' c = J c (wdjaoH a / poJc) 1 ^ 71 • The re- 
sults can be converted to dimensional units by the trans- 
formations J — » JJ' C1 B z — )> B z jioJ r c1 t — » tJ' c /H a and 
E^EpoJ'Jd. 

Shown in Fig. [4] is the result of calculating the distri- 
bution of B z at an applied field of H a = 1 during a field 
ramp at the rate H a = 1 starting from H a = 0. The mag- 
netic flux penetrates from the edges where a critical state 
is established with J « 1. The antidots are clearly seen, 
and the flux distribution reproduces all characteristics of 
the magneto-optical image of Fig. |2j In particular, the 
simulation reproduces the high B z along the array of an- 
tidots, i.e., the flux guidance provided by the patterning. 
The excellent qualitative agreement with the experimen- 
tal image on all visible scales gives strong confidence in 
the correctness of the simulation method. 

Also the d-linei and d-line2, which are clearly seen 
in Fig. [4j are in accordance with both the experiment 
and the Bean model considerations. Interestingly, the 
dark lines d-linei make an angle a = 66° with the strip 
edges, which is slightly larger than expected from the 
Bean model, where Eq. ([!]) gives a = 62° for a/r = 4.4. 
This implies that the current flow is enhanced through 
the bridges between the antidots due to the use of a finite 
creep exponent, n. 

Figure [5] illustrates how the flux distribution evolves 
as the applied field is ramped up. The five panels show 
simulated images from H a = 0.1 to 0.5. Evidently, the 
flux penetration is at all stages greatly advanced along 
the antidot array. In fact, it extends even deeper than 
expected from the Bean model considerations in Fig. Q. 
This is due to the currents flowing in the Meissner state 
part of the film, where they reach the critical value when 
adapting to the constrictions created by the antidots. 

At the field H a = 0.4, a continuous area with non- 
zero B z connects the edges and the central large hole, 
which first at this stage receives sizable amounts of flux. 
Beyond that stage, we find both experimentally and nu- 
merically that one may tune the amount of flux captured 
in the hole by increasing or decreasing H a , and making 
the hole act as a controllable flux reservoir. Note that at 
H a = 0.4 the second large hole is not yet in contact with 
the flux front, but still faintly visible in the figure due 
to its perturbation of the Meissner current flow. In the 
final panel, H a = 0.5, flux is guided from the first large 
hole towards the one below following the connecting ar- 
ray of antidots. Here, the flux motion is parallel to the 
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FIG. 5. (Color online) Distribution of B z in the patterned 
superconducting strip at increasing applied fields of H a = 
O.l(upper), 0.2, 0.3, 0.4 and 0.5 (lower). The inserted color bar 
relates the local magnitude of B z and the image brightness. 



edges. The whole sequence of flux penetration patterns 
during increasing applied field agrees very well with the 
experimental results shown in Fig. [2] 

Figure [6] presents a close-up view of the antidot array 
at H a = 0.4. The upper panel shows the B Z1 where now 
the fishbone structure is seen. The Bean model result, 
Eq. ([2|, predicts for the present array that the d-line2 
forms the angle j3 = 33° with the edge. A dotted line at 
this (3 is included as a guide to the eye in the figure, and 
it demonstrates very good agreement. 
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FIG. 6. (Color online) Distributions of B Z1 J and E along the 
left antidot array, at H a — 0.4. The J-map is superimposed 
with the current stream lines. In the E-m&p, the white spots 
are the bridges between the antidots, where the flux traffic is 
extensive. Color bars relating the image brightness and local 
values are shown. 



The middle panel shows the magnitude of the sheet 
current J with the current stream lines superimposed. 
The main features of the current stream lines are in ex- 
cellent agreement with the Bean model construction of 
Fig. [I] but the curving of the stream lines is less sharp 
due to the finite n. In particular, the cusps creating the 
d-line2 are weak. Note that some current stream lines ex- 
tend into the flux-free Meissner state region, as expected 
in films placed in a perpendicular magnetic field. 

The lower panel shows the magnitude of the electrical 
field E = |E| in the same area. To reveal the overall field 
distribution the images has been plotted on a logarithmic 
scale, with ^-values ranging from 10 ~ 3 to 30. Very large 
fields, E > 30, were found in the bridges of supercon- 
ductor connecting the antidots. Evidently, the antidot 
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FIG. 7. Profiles of (a) B z , (b) J y , (c) \g\, and (d) E y across 
the antidots at H a = 0.4. 



array forms a channel with large traffic of magnetic flux. 
At the same time, the flux motion is much reduced at 
the d-lines, which show up as dark also in the map of 
E, indicating suppressed traffic of magnetic flux. This is 
a generic feature of regions with sharply curved current 
streamlines) 3 -^ 

A quantitative presentation of the profiles for B Z: J, 
g, and E across the antidots array at H a = 0.4 is shown 
in Fig. [7| Clearly, all the quantities are much distorted 
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compared to those of an unpatterned strip pOtHSl j n p ar _ 
ticular, B z has an oscillating behaviour with minima at 
the inner, and maxima at the outer, edges of the antidots, 
where the peak values are comparable to those along the 
strip edge. The sheet current J y is zero within the an- 
tidots and is almost constant in between, with values 
J w 1.1. The the local magnetization g has the shape of 
a step pyramid with flat levels across the antidots. The 
maximum electric field is very high, E > 50, compared 
to a plain strip where E < 1 at the edge, and it is also 
higher than for small indentations at a film edgeP^J 
The high lvalue reflects that all magnetic flux in the 
diamond-shaped region has passed through channels of 
width comparable to the antidot diameter. Numerically, 
the values E w 50 and J ^ 1.1 are consistent, since from 
Eq. ([3), one has 50 1/n « 1.14 for n = 29. 

Previous works concerned with the stability of su- 
perconducting films have found that high electric field 
close to the edge implies that the sample is sus- 
ceptible for ava lanches triggered by thermomagnetic 
instabilities.™! This means that a transverse array of 
antidots, or small defects, are likely nucleation points 
for the instabilities. In thin films, the consequences of 
the instabilities are dramatic, as they often take the 
form of large dendritic structures as observed in many 
materialsPHa 

Finally, note that the results presented in Section II 
were based only on current conservation, and will give a 
quite good description for any sample thickness provided 
the superconductor behaves according to the critical- 
state model. Thus, also the general concept of flux guid- 
ance as presented in this work should be essentially inde- 
pendent of thickness. On the other hand, the simulation 
formalism, which gives a more precise description, was 
derived under the assumption that the thickness is much 
smaller than the size of any lateral structure in the sam- 
ple. Therefore, we expect that the presented simulation 
results will hold only as long as the antidot diameter is 
considerably larger than the film thickness, as was the 
case in the present experiments. 



VI. SUMMARY 

Magnetic flux guidance by linear arrays of antidots in 
type-II superconducting films has been considered exper- 
imentally and theoretically. Experimentally we have used 
MOI to show strong flux guidance in a YBCO film shaped 
as a long strip and patterned with antidots arranged in 
linear arrays. It was also shown that an antidot array 
along the center line can promote flux motion parallel to 
the strip edge. The flux penetration patterns have re- 
vealed that the antidot arrays perturb the overall flow 
of current considerably. In particular, we find that lines 
where the current flow abruptly changes direction, the 
d-lines, agree very well with current stream line patterns 
constructed from the Bean critical-state model. Further 
insight into the flux penetration process was achieved by 
numerical simulations of the electrodynamic response of 
the superconductor subjected to an increasing perpen- 
dicular magnetic field. 

In order to perform simulations of superconducting 
films with complicated antidot patterns it was necessary 
to develop a efficient method for imposing boundary con- 
ditions, both for the sample boundaries and the antidots. 
The simulation of an ascending field ramp produced a 
flux distribution with the same qualitative and quan- 
titative characteristics as found in the experiment. In 
addition, the simulations give a deeper insight into the 
dynamical process of flux guidance by providing informa- 
tion of all electrodynamic quantities at all stages of the 
process. It was shown that the electric field becomes very 
large in a thin channel connecting the antidots, in par- 
ticular close to the edges. Also, maps of the electric field 
show that flux motion is much suppressed elsewhere near 
the antidots. This means that magnetic flux is guided 
into the film via a main route along the antidot array. 
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